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The hydrodynamic behavior of an accretion disk in a close binary system is numerically 
simulated. Calculation is made for a region including the compact star and the gas-supplying 
l- companion. The equation of state is that of an ideal gas characterized by the specific heat 

' ratio 7. Two cases with 7 of 1.01 and 1.2 are studied. Our calculations show that the gas, 

^0 ' flowing from the companion via a Lagrangian LI point towards the accretion disk, forms a 

fine gas beam (LI stream), which penetrates into the disk. No hot spot therefore forms in 
these calculations. Another fact discovered is that the gas rotating with the disk forms, on 
collision with the LI stream, a bow shock wave, which may be called an LI shock. The disk 
becomes hot because the LI shock heats the disk gas in the outer parts of the disk, so that the 
^ , spiral shocks wind loosely even with 7 — 1.01. The LI shock enhances the non-axisymmetry 

QQ ' of the density distribution in the disk, and therefore the angular momentum transfer by the 

' tidal torque works more effectively. The maximum of the effective a becomes ~ 0.3. The 

, 'hot spot' is not formed in our simulations, but our results suggest the formation of the 'hot 

. line', which is the LI shock elongated along the penetrating LI stream. 

o 

§1. Introduction and Summary 

'"q^. 1.1. Discovery of Spiral Structure 

Q ■ In recent years, Steeghs, Harlaftis Sz Horne discovered, with use of the Doppler 

^ , tomography technique, a spiral structure in the accretion disk of the dwarf nova IP 

^ I Pegasi, occurring in its outburst phase. After the discovery of Steeghs et al.,^) 

similar spiral structures have been found on many other accretion disks. ^^"^^ Such 
spiral structures are generally thought to be unobservable during quiescent phases 
^ I of dwarf novae, but Neustroev & Borisov^^ observed them in a quiescent phase in U 

■ Gem. 

Formation of spiral-shape shocks in an accretion disk of a close binary had been 
predicted based on two-dimensional numerical simulations, by Sawada, Matsuda & 
Hachisu, and Sawada et al. ®^ Similar two-dimensional numerical simulations 
performed by many researchers thereafter all confirmed the presence of the spiral 
shock waves. " Self-similar solutions of the spiral shocks were also obtained. 
Angular momentum loss at the spiral shocks provides a different mechanism from 
the conventional a-viscosity one. 

Three-dimensional simulations have also been made on the same object, but 
given no results completely agreeing with one another. Most of the three-dimensional 
numerical simulations of accretion disks so far performed used the particle methods 
such as the SPH method, and mostly failed in finding any spiral shock. 

On the other hand, Yukawa, Matsuda & Boffin observed, with use of the 
SPH method, spiral shocks on the assumption that the specific heat ratio 7 was 1.2. 
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Lanzafame &; Belvedere, Boffin, Haraguchi & Matsuda also observed spiral 

shocks in their recent SPH calculations. This disagreement is considered to be due 
to the difference in resolution, more concretely, in the number of particles. That is, 
use of a sufficient number of particles, which enhances the resolution, leads to the 
observation of spiral shocks. 

With the finite difference/volume method, having a high accuracy compared 
to the SPH method, it is expected that a detailed structure of an accretion disk 
be readily observed. Sawada & Matsuda performed simulation by means of the 
three-dimensional finite volume method for the first time. Although the result shows 
the presence of spiral shocks, no definite statement was made with respect to this 
phenomenon, since the calculation time adopted by them was as short as half an 
orbital period. 

In recent years, our group conducted a series of the three-dimensional finite 
volume calculations and confirmed the presence of spiral shocks in all cases. ' 
These calculations, however, restricted the calculation region to the vicinity of the 
compact star and did not include the companion. The gas was assumed to flow 
into through a rectangular hole located at the Lagrangian LI point. Since most of 
the calculation region was filled with a low-pressure gas, the gas flowing into the 
region from the hole underwent a rapid expansion, thereby forming what is known 
as under-expanded jet. The flow (LI stream) thus formed appeared to interact with 
the accretion disk in a complex manner. The hole, being artificial, made the LI 
stream an undcr-cxpandcd jet having a square cross section, which was unnatural. 
It has therefore been desired that simulation be made with the calculation region 
extending up to the companion, whereby the true shape of the LI stream could be 
observed. 

Bisikalo et al. ' also made very similar works to the present work, and the 
relation between our results and theirs will be discussed in the last section. 

1.2. Difference between 2D and 3D Results 

Although radiative cooling has an important effect on the behavior of an accre- 
tion disk, it is however very difficult to incorporate the effect accurately in three- 
dimensional calculations. To avoid this difficulty, we use as the equation of state of 
the gas the one for an ideal gas and, as the specific heat ratio of the gas 7, a value 
smaller than 5/3, thereby simulating the effect of radiative cooling to some extent. 

What has puzzled us most in interpreting the results obtained by our three- 
dimensional calculations is that they greatly differ from those obtained by two- 
dimensional calculations (see figure 1). With the two-dimensional calculations, we 
can obtain a clear correlation between 7 and the degree of winding-in of spiral shocks: 
smaller pitch angle and tighter winding angle with smaller 7 (Makita et al. , Mat- 
suda et al. ^^)). A smaller 7 corresponds, with the gas undergoing an adiabatic 
change, a smaller temperature elevation. The limit of 7 = 1 corresponds to an 
isothermal change. In the presence of shocks this is not necessarily so, but with 
two-dimensional calculations, a small 7 generally means a low temperature of the 
accretion disk, which corresponds to a low sound speed and thus to a large Mach 
number. It is natural that in this case the spiral shocks wind in tightly. 
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TIME 44.0 STEP 116884 TIME 44.0 STEP 144566 
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Fig. 1. Difference between the 2D and 3D results: The density distributions obtained by 2D cal- 
culation (top) and those by 3D calculation on the rotational plane (bottom), with the specific 
heat ratio 7 being 1.01 (left) and 1.2 (right). The calculation region is restricted to the vicinity 
of the compact star and exclude the companion. With 2D calculations smaller 7 leads to tighter 
winding. On the other hand, 3D calculations do not show such an effect clearly. The strange 
shape of the LI stream is due to the fact that the LI stream comes from a square hole placed 
at the LI point and, as a result, forms an under-expanded jet. (After Makita et al. ^^^) 

On the other hand, as is clear from figure 1, with three-dimensional calculations, 
the obtained spiral shocks seem to wind in very mildly, in particular when 7 = 1.01. 

This finding has long been a mystery to us. Two possible explanations may be 
given for this difference: numerical and physical ones. Numerical explanations are: 

1. The calculation region does not contain the companion, so that the LI stream 
does not have a natural shape. This fact may affect the structure of the accretion 
disk. 

2. In early stages of the calculation, a high-temperature gas is assumed to be 
present over the entire calculation region. This assumption is necessary for 
the calculation to proceed stably. In two-dimensional calculations, an accretion 
disk forms with elapse of time, while the gas initially placed is removed from 
the region of the accretion disk and has no influence on the result. 

3. With three-dimensional calculations, however, the high-temperature gas re- 
mains above the accretion disk and may mix with the disk gas, thereby in- 
creasing the disk temperature. This effect may be significant at the central 
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part of the disk where it is thin. 

In order to solve these problems, the present study enlarges the calculation 
region to an extent covering the companion. Besides, in order to check the influence 
of the high-temperature gas initially placed, the temperature and density of the gas 
are intentionally decreased. As a result it has been found that the above numerical 
factors are not crucial in causing the difference between the calculation results. It is 
then concluded that the difference in the results between two-dimensional and three- 
dimensional calculations should be explained on physical basis rather than numerical 
basis. 

The results of the present calculations taking the companion into account, to 
be shown later, indicate that the basic features of the accretion disk does not differ 
very much from that obtained by calculations without the companion. 

1.3. Penetration of the LI Stream into the Disk and the Loss of the Angular Mo- 
mentum of the Disk Due to the Penetration 

Our calculations covering a region including the companion has revealed that 
the LI stream, on collision with the disk, penetrates into the disk like a spear and 
does not decrease its velocity by forming a hot spot. We refer to this phenomenon 
as "penetration" hereinafter. The presence of the penetration was not convincingly 
confirmed in the previous three-dimensional calculations due to the LI stream hav- 
ing an unnatural shape. However, here, thanks to a visualization technique adopted 
by us, the LI stream reveals its detailed structure. Note that with two-dimensional 
calculations the penetration can never occur, which is understandable from geomet- 
rical consideration. This is an essential difference between three-dimensional and 
two-dimensional calculations. 

The disk gas supersonically rotating around the compact star must form, on 
collision with the high-density, low-temperature LI stream resembling a spear, a 
bow shock, which we call LI shock. The bow shock thus formed may be viewed as a 
hot line or a heated wall rather than a hot spot. At the LI shock the disk gas gives 
the penetrating LI stream its angular momentum. The gas is at last released into 
the inner parts of the disk and mixes with the disk gas. The angular momentum, 
which is transferred from the disk gas to the penetrating LI stream at the LI shock, 
is therefore finally redistributed into the disk, if tidal torque from the secondary 
star dose not work. But, the LI shock and the penetrating LI stream enhance 
the non-axisymmetry in density distribution produced by spiral shocks in the disk 
and, as a result, tidal torque from the secondary star transfers more effectively the 
angular momentum of the disk to the angular momentum of the orbital rotation of 
the binary system. The disk gas therefore loses angular momentum and accretes 
onto the primary star. In figure 2 schematic presentation of the spiral shocks, the 
LI stream and the LI shock are shown. 

Here we present an analogy between a gas centrifuge to enrich uranium and the 
accretion disk with the LI stream. In order to enrich uranium, there have been 
available ultra gas centrifuges, which comprise a cylinder rotating at a high speed 
and containing uranium hexa-fluoride. The cylinder is equipped in the inside thereof 
a stationary tube named "scoop" for recovering enriched or depleted products of 
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Fig. 2. Schematic presentation of tlie flows in the accretion disk, showing: the companion star 
(left end), the accretion disk (a large circle), the spiral shocks, the LI stream flowing from the 
companion to the accretion disk through LI point and the LI shock, on which rotating disk gas 
collides with the LI stream and deflects its direction. 



the uranium compound. The rotating gas forms a bow shock on coUision with the 
scoop, and loses its angular momentum, whereby a meridional flow called counter 
current is formed (Matsuda, Tamura &: Sawada^^^). The scoop works as a sink of 
angular momentum, because the scoop and gas centrifuge are fixed to ground. The 
LI stream in the present case however dose not work as a sink of angular momentum 
itself. The role of the LI stream is only to get angular momentum from the disk gas, 
and the secondary star actually plays the role of a sink of angular momentum. 

When calculation is made by means of a three-dimensional simulation, the an- 
gular momentum of the disk gas is lost on collision with the LI stream and, at the 
same time, the entropy and temperature of the disk gas increase. As a result, the 
sound speed increases and the Mach number in the disk is reduced compared to that 
with two-dimensional calculations, so that the spiral shocks wind in mildly. 

Section 2 briefly describes the method of calculation. Section 3 discusses the 
effect of the numerical assumptions on the results. Section 4 discusses the detailed 
structure of the accretion disk and the LI stream. We compare the three-dimensional 
calculations and two-dimensional ones in detail and shows that the difference can 
be explained on a physical basis rather than on a numerical one. Section 5 gives a 
summary and discussion. 
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§2. Numerical Models 

2.1. Method of Numerical Simulation and Basic Assumptions 

We use, as before, the Simplified Flux Vector Splitting (SFS) finite volume 
method proposed by Jyounouchi et al.,"^^^ Shima & Jyounouchi. ^''^ This method, 
being a variation of AUSM type methods, has the feature of being simple and stable. 
As to the detail of the numerical scheme and the SFS scheme see Makita et al. 
We use a MUSCL-typc technique, thereby keeping the spatial and temporal accuracy 
at second order levels. We consider the g ctS clS cl perfect gas and do not take into 
account such complex effects as viscosity, radiation and magnetic field. The equation 
of state is as follows. 



p = (7 - 1) 



(2-1) 



where p, p, e are the pressure, density and total energy per unit volume of the gas, 
respectively, and 7 is the specific heat ratio; u, v, w are the velocity component in 
the directions of x, y, z, respectively. We use as the specific heat ratio, not the value 
of 5/3, which is that of an adiabatic gas, but smaller values and mainly 7 = 1.01, 1.2. 
In the absence of shocks, the entropy would be kept constant and the temperature 
would change as T oc /0^~^. The temperature change thus becomes nearly isothermal 
at a limit of 7 approaching 1. Our calculations, however, produce shocks, where the 
entropy of the gas increases. Accordingly, 7 ~ 1 does not necessarily mean that the 
accretion disk is nearly isothermal. 

Indeed, with use of the same value for 7, the temperature and the entropy of 
an accretion disk obtained by means of three-dimensional calculations differ from 
those by two-dimensional ones significantly. This is attributable to whether or not 
shocks are present and to the difference in the strengths of the shocks. By selection 
of 7 = 1.01, we intend to increase the cooling effect as much as possible while keeping 
numerical accuracy. The larger 7 = 1.2 is selected because a 7 larger than this would 
not produce an accretion disk itself. We thus study the two extreme cases. 

We consider, as the model of a binary, a system consisting of a compact star 
(primary star) with mass Mi and a mass-losing companion with mass M2 ■ The mass 
ratio q = M2/M1 is fixed at 1 in the present study, while with IP Peg q ~ 0.5. As to 
other mass ratios see the paper by Matsuda et al.,^^) although they restricted their 
computational region to the vicinity of the primary star. 

All physical quantities are normalized. The length is scaled by the separation A. 
With the rotational frequency of the binary being H, the time is scaled as l/f2, so 
that the orbital period becomes 2-it. The density at the surface of the companion is 
taken to be 1. The gravity constant G is eliminated through the following equation. 

= G(Mi +M2). (2-2) 

A Cartesian coordinate system with the origin located at the center of the pri- 
mary star is used. The x-axis is a line combining the centers of the primary and the 
companion. The center of the companion is at (—1, 0, 0). The z-axis is perpendicular 



Numerical Simulation of an Accretion Disk. 



7 



Table I. Parameters. 
Model Co Pi 5 ci 7 



A 


0.02 


10" 






1.01 


B 


0.02 


10" 


-5 


Vw 


1.2 


C 


0.02 


10" 


-7 


0.02 


1.01 


D 


0.02 


10" 


-7 


0.02 


1.2 


E 


0.1 


10" 


-5 


Vw 


1.01 


F 


0.1 


10" 


-5 


Vw 


1.2 



to the orbital plane and is oriented in the same direction as the angular momentum 
vector of orbital rotation. The x — y plane therefore constitutes the orbital plane. 

The calculation region is a rectangular parallelepiped covering from (—1.5, —0.5, 0) 
to (0.5, 0.5, 0.5). With the assumption of symmetry of physical quantities in relation 
to the orbital plane, the calculations are made only in the region of upper half. This 
region is divided with grid points of 201 x 101 x 51. 

2.2. Boundary Conditions and Initial Conditions 

The primary star is assumed to be a hole occupying a region covering 3x3x2 
and centered at the origin. The inside of the primary star is always filled with a 
gas having a density, sound speed and velocity of pi, c\ and 0, respectively. The 
pressure in the inside of the hole is kept sufficiently low in order to always absorb 
the surrounding gas. Here p\ and c\ are parameters. The mass accretion rate onto 
the main star is obtained by solving Riemann problem at each time step. 

The companion is assumed to fill the critical Roche lobe. The inside of the 
companion is filled with a gas with density po = 1) sound speed cq and velocity 0, 
while the gravity is neglected in the companion. Here co is taken as a parameter 
and its effect is studied. The assumption that the gas inside the companion has 
velocity does not mean that there is no gas outflow from the companion. The gas 
does flow out when the pressure outside the companion is lower than that inside the 
companion. The outflow flux is calculated by solving Riemann problem at each time 
step. 

The outside of the outer boundary is assumed to be always filled with a gas 
having a density pi, sound speed c\ and velocity 0. Use of this boundary condition 
insures stable calculation. Here also inflow or outflow of the gas through the outer 
boundary is possible and the flux is calculated by solving Riemann problem between 
the inside and outside of the calculation region at each time step. 

As the initial conditions, the whole region except the inside of the companion is 
occupied at time t = by the gas with density pi, sound speed c\ and velocity 0. 
The initial density p\ and initial sound speed c\ are parameters, the effects of which 
are studied later. The model parameters used in the present work are summarized 
in the table 1. 

Typically, calculations are continued up to an orbital period of more than 10. 
We confirm that the disk mass does not change appreciably after a few rotation 
period, so we may say that our computational results are nearly in steady state 
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except insignificant oscillations. 

§3. Effect of Various Parameters on the Results of Calculations 

3.1. Effect of the Conditions of Gas Placed Initially in the Calculation Region 

In our previous calculations, " a gas with a density pi = 10~^ and a sound 
speed ci = -\/lO was placed at time t = 0. This high sound speed means that the 
gas has a high temperature. This selection was made, since it was considered that 
such a high-tcmpcrature, low-density gas hardly falls onto the accretion disk or, if 
falls at all, affects the accretion disk only to a small extent. 

The results obtained by the three-dimensional calculations, however, differed 
greatly, as described in Introduction, from those by the two-dimensional calcula- 
tions. That is, the spiral arms with the three-dimensional calculations wind in more 
mildly than those with the two-dimensional calculations, even with the specific heat 
ratio 7 = 1.01. This is attributable to the high temperature of the accretion disk 
obtained in the three-dimensional calculations, which in turn may be due to the 
high-temperature gas initially placed. 

Unlike with two-dimensional calculations, the gas initially placed will, during 
progress of calculation and after an accretion disk has been formed, remain above the 
accretion disk. This initial gas may mix with the accretion disk, thereby increasing 
the temperature. We selected the density of the initial gas to be 10~^, which we 
considered to be sufficiently small. However, the central part of the accretion disk, 
having a small thickness, may have been affected. 




Fig. 3. Effect of the initially placed gas: Iso-density surfaces at log p = —3.3 of the disk for Model A 

with high temperature initial gas (a) and that at log p — —4.0 for Model C with low temperature 
and low density initial gas (b). The Model C shows clearer spiral shocks. In both cases, spiral 
shocks wind more loosely than that in two-dimensional calculations with 7 = 1.01. 

In the present study, we compare two cases with the initial gas having different 
densities and temperatures, in order to check their influence. Figure 3 shows the 
results. We may observe that Model C with low temperature and low density initial 
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gas shows clearer and more openly wind spiral shocks than Model A with high 
temperature and high density initial gas. However, the difference may not be very 
significant, and in both cases, spiral shocks wind loosely even with 7 = 1.01. 

The key role is played, not by the temperature itself, but by the heat capacity, of 
the initial gas. A sufficiently low density of the initial gas therefore causes no problem 
with respect to heat capacity even when it has high temperature. Therefore, we may 
conclude that our previous choice of the initial condition is good enough, as far as 
the cases of 7 = 1.01 are concerned. The cases of 7 = 1.2 are slightly affected by the 
initial condition as will be discussed in 3.4. 

3.2. Effect of the Surface Temperature of Companion 

In our previous calculations, the gas placed at LI point was assumed to have a 
sound speed of cq = 0.1. The present study includes the case of the gas present on 
the companion having a sound speed of Cq = 0.1, which is, however, very large in 
actuality. If the binary has a typical speed of 500 km s^^, the sound speed becomes 
50 km s~^, corresponding to a temperature of as high as 2.5 x 10^ K. 

This high sound speed has been selected, because a high temperature on the sur- 
face of the companion increases the rate of the gas flowing out through the LI point, 
so that calculation can proceed fast. Since, however, the above high temperature is 
unrealistic, the present study also includes the case of cq = 0.02 corresponding to 
10^ K. 

Figure 4 shows, with 7 = 1.01, the influence of the surface temperature of the 
companion. With the high temperature, the LI stream becomes thick and its inflow 
rate increases. On the other hand, with the low temperature, which may be more 
realistic, the LI stream becomes thin. Model E shows clearer and more widely open 
spiral shocks than Model A, because of larger disk size. 



(a) (b) 




Fig. 4. Effect of the surface temperature of companion: Iso-density surfaces at log p = —2.5 for 
Model A with high surface temperature of companion (a), and that at logp = —3.3 for Model 
E with low surface temperature (b). With the high surface temperature of companion, the LI 
stream becomes thick and the inflow rate high and the disk becomes large. In Model E, spiral 
shocks become clearer and wind more loosely than that in Model A. 



10 



H. Fujiwara, M. Makita, T. Nagae and T. Matsuda 



3.3. Effect of Specific Heat Ratio 

As commented in Introduction, the influence of radiative cooling is very impor- 
tant in accretion disks. Without this cooling, the temperature of an accretion disk 
would become the virial temperature, which is very high. It is however not practical 
to perform any three-dimensional simulation of radiative fluid dynamics accurately 
with computing capability currently available. 

Therefore, workers in this field have been used either an eqation of state with low 
specific heat ratio like us, that of iso-thermal gas or that of iso-entropy gas. A low 
specific heat ratio is selected in the present work in order to emulate the influence 
of radiative cooling. Although the ionized plasma gas present in actual accretion 
disks has a specific heat ratio of 5/3, the present calculations use 7 = 1.01, 1.2. 
Our calculations thus, in a strict sense, simulate a poly-atomic molecular gas. In 
this case, when the gas is compressed, its energy, being distributed to the degree 
of internal freedom, will not contribute to translational motion to a large extent, so 
that the temperature elevation is suppressed. This internal energy is comvected with 
gas and eventually accreted by the central object, and is not radiated away to the 
outer space. Nevertheless the low ratio of specific heats emurates radiative cooling 
somehow. The caluculation including the radiative transfer is now work in progress 
(Nagae et al. 2001). 

In the limit of 7 = 1, the gas becomes isothermal. The equation of state adopted 
by us, however, does not allow 7 = 1. We therefore select 7 = 1.01, which is the 
minimum value in view of calculation accuracy. As a counterpart, we select 7 = 1.2, 
above which experience tells us no accretion disk can ever form. Figure 5 shows the 
results. 

With 7 = 1.01, a clear accretion disk and shock waves appear under all the 
initial conditions or boundary conditions assumed. On the other hand, with 7 = 
1.2, the shape of the accretion disk finally obtained varies depending on the initial 
conditions and boundary conditions. In Model B and F the accretion disks are highly 
distorted, while in Model D we may observe a normal accretion disk. In all cases, 
the temperature of the disk is more higher than that of the disk in its outburst. 

Bisikalo et al. ^^^"^^^ argue that no accretion disk forms with 7 = 1.2, which 
is not completely inconsistent with our results. At any rate, it is expected that a 
specific heat ratio larger than 1.2 may produce no accretion disk under any initial 
or boundary conditions. 

§4. Detailed Structure of Flow 

4.1. Penetration of the LI Stream into the Accretion Disk and Formation of an LI 
Shock 

In order to investigate the detailed structure of the accretion disk and the LI 
stream, let us investigate the details of some specific models in this section. Since 
the case of specific heat ratio 7 = 1.2 corresponds to too high a temperature of the 
disk, only the case 7 = 1.01 is studied here. In the present case, cq = 0.02 is taken 
as the sound speed on the surface of the companion. 
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(a) W fcj 




Fig. 5. Iso-density surfaces for a specific heat ratio of 7 = 1.2, under different boundary conditions 
and initial conditions: (a) Model B at logp = —3.7. (b) Model D at logp = —4.0, with low 

density and low temperature of initial gas. (c) Model F at logp = —2.9, with high temperature 
at the surface of the companion. With 7 = 1.2, disk is affected significantly by both initial and 
boundary conditions. 



Figure 6 shows the iso-density surfaces at various levels of Model C. And figure 
7 shows the iso-density surface of the same model from a different view angle. One 

may observe cross-sectional structures of a spiral shock and the LI shock. The spiral 
shock seems to be concave to the upstream, while the LI shock is convex. One may 
conclude that the LI shock is a bow shock. 

The figures so far presented show only iso-density surfaces, which cannot reveal 
any detailed structure, in particular the velocity distribution, of the flow. Hereafter 
let us investigate the interaction between the LI stream and the accretion disk, with 
use of the velocity distribution. 

Figure 8 shows velocity distributions on x-y planes at various heights by means 
of the Line Integration Convolution (LIC) method. The figure shows that on the 
orbital plane of 2; = 0, gas in the accretion disk collides with the LI stream and, 
as a result, loses its angular momentum and falls towards the primary star. A bow 
shock is then formed close to and parallel with the high-density, spear-likc LI stream 
(above the LI stream in the figure). On z = 0.08, the gas of the accretion disk rotates 
about the central star. The influence of the LI stream is thus limited to a region 
near the orbital plane. 

Figure 9 presents, in order to observe the above phenomena more closely, flow 
patterns, drawn by means of the LIC method, on the iso-density surface. The figure 
shows the internal structure more clearly with increasing density of the iso-density 
surface. Figure 10 shows a three-dimensional view of the streamlines. 

These pictures are particularly interesting, because they show two remarkable 
features. 1) The LI stream penetrates into the accretion disk. 2) The disk gas having 
collided with the LI stream abruptly changes, due to the LI bow shock formed, its 
direction towards the primary star. 
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4.2. Influence of the LI Stream on the Disk Gas 

Figure 11 shows, in order to study the influence of the LI stream on the disk 
gas, the distribution of the entropy related function of the gas S = p/ obtained by 
two-dimensional and three-dimensional calculations, with the abscissa representing 
the distance from the central star and the ordinate the entropy: log(S'). Dispersion 
in the plots in the ordinate direction corresponds to that in the values at various 
angular positions. 

We discuss, in particular, the case of 7 = 1.01, where, S being almost the 
same as temperature, the figure can be considered to be showing the temperature 
distribution. The figure clearly shows a smaller dispersion with two-dimensional 
calculations than with three-dimensional calculations. This means that the entropy 
and temperature distributions in two-dimensional cases are nearly axisymmetric, 
which in turn means that the entropy increases due to spiral shocks, only to a 
moderate extent. The entropy values in the central zone (the left end) and the 
peripheral zone (the right end) arc lower with two-dimensional calculations than 
with three-dimensional ones. A large dispersion in the values with three-dimensional 
calculations is caused by the LI shock generated on collision between the LI stream 
and the disk gas. In particular, a horizontal linear distribution of S, seen in a 
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Fig. 7. The same as figure 6 at log p = —4.0 but seen from a different angle. One may observe cross 
sections of a spiral shock and the LI shock. The spiral shock is concave to the upstream, while 
the LI shock is convex to the upstream. This means the LI stream is a bow shock. 



region corresponding to large radii, is due to the LI shock. Accordingly, with three- 
dimensional calculations, a high entropy or high temperature of the gas decreases 
the Mach number of the gas, so that spiral shocks wind in mildly even with 7 = 1.01. 

Figure 11 also shows the entropy distribution for the case of 7 = 1.2. In this 
case, although the distribution is wider in three-dimensional calculations than two- 
dimensional ones, the absolute values are almost the same in both calculations. The 
angle of winding-in of spiral shocks is therefore almost the same with both two- 
dimensional and three-dimensional calculations. 

4.3. Effective a 

The accretion disk has a large non-axisymmetry in its density distribution be- 
cause of the loosely wound spiral shocks in our calculations. In the spiral shock 
model of angular momentum transfer, tidal torque from the secondary star transfers 
the angular momentum of gas in the accretion disk to the angular momentum of 
the orbital motion of the binary system. Angular momentum transfer becomes more 
effective as spiral shocks wind more loosely and, as a result, the non-axisymmetry 
of density distribution becomes large. 

In our calculations, the disk gas loses its angular momentum not only through 
spiral shocks, but on interaction with the penetrating LI stream. 

The disk gas loses its angular momentum on collision with the penetrating LI 
stream. In this process, angular momentum is transferred directly from the disk 
gas to the penetrating LI stream. This asngular momentum transfer is extremely 
effective on the equatorial plane. In figure 8, the streamlines of the disk gas is bent 
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Fig. 8. Flow patterns on the orbital plane and on x-y planes at various heights, for Model A, are 
shown by means of the LIC method. The degree of whiteness represents the magnitude of 
density. From top to bottom z — 0, z — 0.04; and z = 0.08. On the orbital plane with z = 0, 
the gas in the accretion disk collides with the LI stream and, as a result, loses its angular 
momentum and falls towards the central star. On z = 0.08, the gas of the accretion disk rotates 
about the central star. The influence of the LI stream is thus limited to a region near the orbital 
plane. 



on collision with the penetrating LI stream. This mechanism is effective, but is 
closed in the disk. The angular momentum transferred to the penetrating LI stream 
is finally released into the central part of the disk. This mechanism is, therefore, 
angular momentum re-distribution rather than angular momentum transfer. 

The LI shock produced on collision between the disk gas and the penetrating LI 
stream contributes to a non-axisymmetry in density distribution in the disk, and the 
disk gas loses its angular momentum through the LI shock. These two mechanisms 
are effective in the peripheral part of the disk into which the LI stream is penetrating. 
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Fig. 9. Isodensity surfaces and flow patterns thereon, together with a flow pattern on the orbital 
plane, of Model A drawn by LIC method. Prom top to bottom at logp = —3.3, —3.0, —2.7 all 
at time t = 75.36. 



We calculate the effective a, according to Frank, King & Raine, 



M 



a = 



(4-1) 
(4-2) 



where v is kinematic viscosity, M is mass accretion rate, S is surface density, i?* is 
radius of the primary star, R is distance from the primary star, Cg is sound speed 



and H is scale height. We neglect ^R^/Rj because R*/R <^ 1. Figure 12 shows 

the effective a at each orbital period for Model A. The effective a varies with both 
of radius and time, its maximum being a ~ 0.3. 

In Spruit and Livio &; Spruit, the effective a depends on the mass ratio 
and the pitch angle of spiral shock. However, the analytical value of the effective a 
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Fig. 10. Three-dimensional view of streamlines of model A. 



is of order of ~ 10~^, while in two-dimensional calculations by Matsuda et al. -"^^^ 
the effective a is of the order 0.1 with large 7. 

The effective a is large in the peripheral part. In this region, the non-axisymmetric 
density distribution is enhanced by the LI shock, and therefore, the disk gas loses 
effectively its angula momentum due to tidal torque from the secondary star. The 
disk gas also loses its angular momentum on collision with the penetrating LI stream 
in the peripheral part. On the other hand, in the central part, the effective a is of 
the order 0.1, and this value is similar to that in two dimensional calculations with 
large 7 and loosely wound spiral shocks. 

4.4. Formation of the 'hot spot' 

In our three-dimensional simulation, the effect of radiative cooling is neglected, 
and therefore the accretion disk becomes hot. This hot disk becomes thick and of 
low density, so that the dense LI stream can penetrate into the disk and the 'hot 
line' is formed along the penetrating LI stream. 

On the other hand, in two-dimensional simulations, the disk cannot expand along 
the z-axis, and therefore the density of the disk becomes larger than that of the LI 
stream. And the LI stream collides with the disk edge and does not penetrate into 
the disk. "^^^ Rozyczka & Schwarzenberg-Czerny and Rozyczka^^^ performed 
two-dimensional simulations of the collision between the LI stream and the disk 
edge. They showed that the collision between the LI stream and the disk edge forms 
two shock waves and concluded that these shock waves are considered as the 'hot 
spot'. 
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Fig. 11. Radial distribution of entropy log(p/p''') in the accretion disk. The top row corresponds 
to 7 = 1.01, and the bottom row to 7 = 1.2. The left-side column shows the case with two- 
dimensional calculations and the right-side column that with three-dimensional calculations. 
Other parameters are pi = 10~'^,ci = 0.02, and co = 0.02. In three-dimensional calculations, 
entropy increases rapidly at the peripheral zone due to the LI shock. This increase of entropy 
is significant with 7 = 1.01. 



The difference of tlie density between tlie disk and the LI stream is important. 
Whether or not the LI stream penetrates into the disk depends on the magnitude of 
each of the densities. If the density of the disk is larger than that of the LI stream 
in the three-dimensional simulation, the LI stream will be stopped on collision with 
the disk edge and the hot spot will be formed. 

In order to investigate collision between the LI stream and the dense accretion 
disk, we suddenly changed the density of the surface of the secondary star from 1 to 
0.5, after the accretion disk has reached steady state. In this case, the density of the 
LI stream decreases and becomes of the same order of magnitude with the density 
at the disk edge. Figure 13 shows the density contours on the equatorial plane in the 
vicinity of the collision point. The LI stream is stopped on collision with the disk 
edge, and shock waves arc formed at the collision point. These features are similar 
to the results of two-dimensional simulations. 

If the effect of radiative cooling is considered, the accretion disk will becomes 
cool, thin and of high density. In such conditions, we consider that the LI stream 
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cannot penetrate into the dense disk and that the 'hot spot' is formed on cohision 
between the LI stream and the disk. 

4.5. Flows on the Companion 

We have so far paid attention to the accretion disk and the LI stream. Since 
our calculation region includes the companion, it is also of interest how the flows 
on the surface of the companion behave, although we may not recognize this with 
observational means currently available. 

The coordinate system we employed is Cartesian, which can hardly represent 
precisely the companion having a nearly spherical surface. Outer boundaries may 
also affect the surface flow of the companion. We have to examine the results with 
these restrictions in mind. Figure 14 shows flows on the surface of the compan- 
ion for case A at 9 orbital periods. Figures are produced using Pixcel Exposure 
Method (PEM). Worth noticing is that flows towards the LI point start even at 
high-latitude regions. This main flow pattern is almost steady, although small flow 
patterns change. 

It is also interesting to note that eddy patterns have been formed on the surface 
of the companion. This is considered to be similar to vortices of low pressure or 
high pressure generating on the surface of Earth or Jupiter. It is expected that flows 
on the surface of the companion in a binary system do not form any simple regular 
pattern but become complex flows. However, the above restrictions allow no definite 
conclusions to be drawn. Calculations with generalized curvilinear coordinates fitting 
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Fig. 13. Density contours on the equatorial plane. The density of the LI stream at the collision 
point is of the same order of magnitude with the density of the disk edge. The LI stream is 
stopped on collision with the disk edge and shock waves are formed. These features are similar 
to the results of two-dimensional simulations. 

the shape of the companion surface are desired. 

§5. Conclusions and Discussion 

5.1. Conclusions 

We perfomed numerical simulations for an ideal gas constituting an accretion 
disk in a close binary with mass ratio 1, with a specific heat ratio of 7 = 1.01, 1.2 
and compared two-dimensional with three-dimensional calculations. Prom the results 
obtained, we conclude as follows. 

1. Whether or not the calculation region may include the companion, no significant 
difference appears in the results. 

2. An initial condition for gas to be placed at t = of 7 = 1.01 does not have a 
large influence, while that of 7 = 1.2 has. 

3. A high temperature of gas present on the surface of the companion gives a thick 
LI stream, while a low temperature a fine one; both do not have very much 
influence on the structure of the accretion disk. 

4. A specific heat ratio of 7 = 1.01 gives, for any parameter, a neat-shaped ac- 
cretion disk with spiral shocks. On the other hand, with 7 = 1.2 the accretion 
disk formed sometimes takes a greatly distorted shape depending on the values 
of parameters. 
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Fig. 14. Flows on the surface of companion produced by Pixcel Exposure Method for case A at 
t = 56.52. The value of density on which the flows are plotted is logp = —3.5. Gas which 
outflows from high-latitude regions flows toward the LI point. This main flow pattern is almost 
steady, although small flow patterns change. Much attention should not be paid on vortices 
at the corner of the rotational plane, because they may be affected by the outer boundary 
condition. 



5. The LI steam does not reduce its velocity on collision with the accretion disk 
or form a hot spot, but forms a bow shock, which can be called an LI shock, 
or a hot line rather than a hot spot. 

6. The LI shock and the penetrating LI stream enhance the non-axisymmetry in 
density distribution in the disk, so that angular momentum transfer by tidal 
torque from the secondary star becomes more effective. 

7. Complex vortices may be present on the surface of the companion. 

5.2. Discussion 

5.2. L Bisikalo et al.'s work 

Attention should be paid here to a series of calculations made by Bisikalo et 
g^j_32)-36) 'pj^gy performed, like ourselves, numerical simulations of accretion disks 
by means of the finite difference method. Our calculation scheme differs from theirs in 
that we use a far larger number of cells, thus achieving a higher calculation accuracy, 
although our calculation region is a little narrower than theirs. Nevertheless, their 
results, varying to some extent though, are basically similar to ours. In particular, 
their results, like ours, show no hot spots. They further argue that the luminosity 
curve of a close binary can be explained even in the absence of hot spots. 
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We also differ from them in the way to explain the mechanism involved in gener- 
ation of shocks. We contend, as before, that spiral shocks are formed due to the tidal 
force exerted by the companion, and that the LI shock occurs on collision between 
the LI stream and the disk flow. On the other hand, they do not admit the effect of 
the tidal force and, further, propose a mechanism of the shock formation as caused 
by collision between the disk flow and the envelope. The fact that our results are 
similar to theirs is still worth noticing, although the interpretations are somewhat 
different, and desired to be checked by a third party. 

5.2.2. Bath et al.'s work 

Bath, Edwards k. Mantle ^^-^ have discussed the penetration of the LI stream 
into the accretion disk. According to their discussion, whether or not the LI stream 
penetrates into the accretion disk depends on the magnitude of each of the densi- 
ties. That is, an LI stream having a density larger than that of the accretion disk 
penetrates into the accretion disk. Our calculations just correspond to this case; the 
density of our LI stream is as high as about 10 times that of the accretion disk. On 
the other hand, an LI stream having a low density will form a hot spot on collision 
with the accretion disk. This case probably corresponds to the quiescent phase hav- 
ing a hot spot. In this case the disk temperature is lower than that in the outburst 
phase due to the lack of the LI shock, and the spiral shocks wind more tightly and 
they may not be observable by the Doppler tomography technique. 

Bath ct al.^'''^ call the LI shock, or the hot line wc refer to in the present paper, 
which we have found by means of calculations, a shock-heated wall. With respect to 
observation of an actual penetration, Chochol et al. argue that it is observable in 
the symbiotic star CI Cyg, and Skidmore et al. also argue in the dwarf nova WZ 
Sge. 

The magnitude of the density of an LI stream depends on the surface temper- 
ature of the companion. With a low surface temperature, the LI stream becomes 
thin and of low density, so that it may not penetrate into the accretion disk. 
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